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We discuss the orbital magnetism and the Hall effect in the weak magnetic field in two dimensional 
Dirac fermion systems with energy gap. This model is related to the graphene sheet, organic 
conductors, and d-density wave superconductors. We found the strong diamagnetism and finite 
Hall conductivity even in gapped systems. We also discuss the relation between the weak-magnetic 
field formalism and the Landau quantization with the Euler-Maclaurin formula. 
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Introduction — Experimental studies of graphene have 
revealed interesting physical propertiespllsj. An anoma- 
lous quantum Hall effect is observed where the quan- 
tized conductance is described by a xy = ±4j-(4n + 2) 
(n = 0, f , 2, • • • ). This result is explained by the band 
structure of graphene which has gapless point due to the 
hexagonal lattice and a linear dispersion around the gap- 
less point d, 0, Moreover, the minimal conductivity 

^min — 

e 2 /h par channel is observed, but theoretical ex- 
planation for this is not successful 0,1a, 0,01- If the Fermi 
energy is near the gapless point, the low-energy property 
of this system is described by the two-dimensional (2D) 
Dirac fermions. Similar gapless system is also found in 
the organic conductor cv-(BEDT-TTF)2l3, which has a 
tilted linear dispersion [1, @|. Moreover, the high-T c su- 
perconductor with ri-density-wave gap is also discussed 
as a Dirac fermion system 13, U [ . 



The physical properties of Dirac fermions in magnetic 
fields have been studied by assuming the Landau quan- 
tization. On the other hand, quite recently, Fukuyama 
discussed the magnetic susceptibility and the Hall effect 
of the gapless Dirac system using the formalism of the 
weak magnetic field [12J. It is predicted that this sys- 
tem shows a strong diamagnetism. On the other hand, 
gapped systems are also interesting. Opening of a en- 
ergy gap is also studied experimentally. The physical 
meaning of gap is directly related to the electron density 
imbalance between the two sublattices of the bipartite 
hexagonal lattice of graphene. We are also interested 
in the relation between the results of the theory of the 
strong magnetic field where the Landau quantization is 
essential and the weak magnetic field, especially for the 
magnetic susceptibility. In this paper, we generalize the 
Fukuyama's work[T2] to the gapped systems, and discuss 
the relation to the strong magnetic field. 

The single particle Hamiltonian of the 2D Dirac 
fermion is given by 



Ho = v(p x a x + PyO-y) + a z A. 



(1) 



where v is the Fermi velocity, A is the energy gap, and 
<7j (i = x, y, z) is the Pauli matrix. In the presence of the 



magnetic field, the momentum operator p = -lfiV in the 
Hamiltonian is replaced by II = p+eA/c where V x A = 
(0,0, B). The Hamiltonian in the many body system is 
given by the field operator ^(r) = [ip\(r),^_(r)], 



H= d 2 r¥(r)H t>(r) 



(2) 



Orbital Magnetism — The general formula of the 
magnetic susceptibility of the interband system in 
terms of the temperature Green function is derived by 
Fukuyma[l3j]. This result is obtained by the Luttinger- 
Kohn representation 14) for the basic functions and 
the Fourier expansion of the vector potential A(r) = 
-iA q {e iq r - e- iq r ) with B = q x A y 



QyAq as, 



(3Vc 2 h 2 



(3) 



n k 



where (3 = and V are the inverse temperature 

and the volume of the system, respectively. The tem- 
perature Green function is given by G = G(k, iu) n ) — 
(iu) n - Ho/h)^ 1 , where iu) n = iuj n + \p + mga(u}„)T]/h, 
with oj n = (2n + \)ir/ fih being the Matsubara frequency 
of fermions. Here, we have introduced the scattering rate, 
r = — Im E R neglecting the frequency dependence, for 
simplicity. The matrix trace tr [■ • • ] of eq. is calcu- 
lated as 



16v 4 k 2 k 2 y 



[(iuj n ) 2 - v 2 k 2 - A 2 ] 4 [(id>„) 2 - v 2 k 2 - A 2 } 2 



(4) 



where A = A/h. Integrating out the wave number fc, we 
have 



Tr En- ©] 



* " J 2,TTV 2 [{lij n ) 2 - A 2 ] 



= F x {iQ n ). (5) 



Finally, taking the cj„-summation as the contour integra- 



2 



eV 1 



tion, the susceptibility is obtained as 

'Fx(iu> n ) (6) 
duf(hw) + if) - - if )] , 



2 4 

e u 



1 



c 2 ft 2tti 



(7) 



where /(a;) = (e^( x_ W + is the Fermi distribu- 

tion function and V = r/7i. Especially, in the zero- 
temperature limit f(x) — > 6([i — x), we have 



6ir 2 c 2 h 2 A 

2 2 

e v 



Tan 

r 



37r 2 C 2 7l 2 T 2 + /i 2 



r 

(A 



Tan 



0). 



(8) 
(9) 



In the A — > limit, this result coincides with that of the 
gapless systemjl^], as expected. Moreover in the clean 
limit r — > 0, this gives the 5-function with negative sign. 
The chemical potential dependence of the susceptibility 
is shown in Fig. QJa). This result shows that a strong 
diamagnetism appears even in the presence of a gap, and 
the susceptibility takes the minimum value when /i is in 
the middle of the gap. This result is similar to the case 
of bismuth [l3j|. 

Now we will show that the magnetic susceptibility 
is also obtained by the Landau quantization formalism. 
To do this, it is sufficient to derive eq. ([5|). The Hamil- 
tonian (TT]) in the magnetic field is given by 



Ha 



A vIL- 
vIU -A 



(10) 



where II -t = H x ± illy . Since the commutation relation 
between these operators is [n±,n+] = ^f2eBH/c, there 
is the correspondence to the creation and the annihila- 
tion operators of the harmonic oscillator: 11+ — > ^/2ja\ 
II_ -> V2ja where I = yJch/eB and eB > 0. Then the 
eigenvalue and the eigenstate of (fTT)|) is obtained as 



Hl\k))=M 2 k \k)), M k 



h 2 v 2 



I 2 



(fc>0), |0» = 





|0) 



(11) 

(12) 



where k = 0, 1,2,3, ••• denote the Landau levels with 
positive energy, and \k) corresponds to the number state 
of the harmonic oscillator. Note that k = state is spe- 
cial: one component has zero amplitude, and the normal- 
ization factor is different from k > cases. 

According to the functional integral method, the ther- 
modynamic potential is given by the temperature Green 



function as 



^ oo 

n(B) = - — J2 Tr l °s(-^n + Ho/h) 



2V 



n— — oo 

oo 



n— — oo 



J2 53 logics) 2 - (M k /ny 



k=l 



+ -log|(iw„) 2 -(A/ A) 2 | 



(13) 



Here factor V/irl 2 in eq. (| 13|) is the degeneracy of the 
Landau levels for |fe| > 0. For k — 0, the degeneracy is 
V/2tt1 2 . Using the Euler-Maclaurin formula, 



1 — r°° i 

- 5 (0) + ^#)^ / g(x)dx--g'(0) ) 

1 k=i J ° 11 



(14) 



the summation over the index of the Landau levels in 
eq. (|13|) is evaluated in the B — > limit as 



1 

ri 2 



1 



dx log 



(iu> n )' z - 2xv 2 /l 2 - A 



Girl 4 (iCj n ) 2 - A 2 



(15) 



Since the first term of eq. (|15|) is constant with respect 
to the magnetic field, the contribution to the magnetic 
susceptibility is given as 



1 / ch\ 2 d 2 n 



2V V ev 2 J OB 2 3ttv 



1 OO 

— V 



1 



n— — oo 



{iw n ) 2 - A 2 



(16) 



Thus eq. (O is obtained. This means that the two differ- 
ent formalisms for the magnetic susceptibilities give the 
same result. 

Hall effect — Next, we consider the transport proper- 
ties of this system. The conductivity of the 2D Dirac 
system is already derived in Ref. llj. To clarify the 
derivation of the Hall conductivity a xy , we rederive this 
result first. The conductivity is given by the Kubo for- 
mula, 



Re a.,,. — lim 



Im n^(0,o; + it?) 



(17) 



where n^(q,w) = n^(q,w) - IL_„,(g, 0). The polariza- 
tion function in the Matsubara form is given by 



1 ( (in 

lW<jr,iO = - / dre 1 ^ (^^(9,^^(0,0)), 
v Jo 



11^(0, ii/ m ) 
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(18) 
(19) 
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FIG. 1: (a) Magnetic susceptibility % = x3i" 2 c 2 ft 2 /e 2 v 2 , 
(b) Conductivity a = o2-Kh/e 2 , and (c) Hall conductivity 
a xy = a xy An 2 c/e 3 v 2 B of the 2D Dirac fermions in the weak- 
magnetic field, as functions of the scaled chemical potential 
(x/T. The results (a), (c) of gapless case (A = 0) are derived 
in Ref. [H. 



where G+ = Q(k,Lu n + v m ) with v m = 2nm//3H being the 
Matsubara frequency of bosons. The current operator is 
given by 

J^(q,r) = -evY / ^ t (k + ,T)a^(k-,r), (20) 



where = k±q/2. The matrix trace tr [■ ■ ■ ] in eq. |T9 
for /i = v — x is calculated as 



2[(kD„)(io}+) - A 



[(id;,,) 2 - v 2 k 2 - A 2 ][(id>+) 2 - v 2 k 2 - A 2 ] 



F 2 (iL) n ,iQ+), 



where ujj~ = lu„ + v m . The summation of the Matsubara 



frequency is carried out in the following way 

du>f(huj) 



i V Fi(iu> n , iu>+) = — !- 

pfl '-^ 27T1 



(22) 



x [Fi(uj + iT,uj + ii> m 
+ Fi(uj — \v m — if, w 



-iT) 
if) 



-Fi(u 
F l {u- 



lT, LO 



iF, ui - 



-J) 
if)]. 

After expanding eq. (|22[) in the linear order of iv m and 
the partial integration of w, it follows from eqs. (JT7J) and 
(Un]) that the conductivity is calculated as 

" = 4n / <lj\-f'(hw)}A 2 (hij,T,A), 



4tt 2 

a 2 (x,t,a) 
r 2 - a 2 + 

2fbj 



(23) 



Tan 



_rA 2 +r 2 



2rb 



+ 1. (24) 



Here the last term in eq. (f2~4"|) stems from the first and 
forth terms in eq. (|22[) where two frequencies have the 
same analytic properties, while the other terms originate 
from the second and the third terms. At zero temper- 
ature, f'(x) = —5{x — /z), we have tr = t^Aq([j,,T, A.) 
which is derived in Ref. [llj. The chemical potential de- 
pendence of the conductivity is shown in Fig. [ljb). In 
the limit fj,/T — ► 0, we have the minimal conductivity, 
0min = + A 2 /r 2 ) -1 . The dependence of the mag- 
netic field B can be calculated in the similar way of a xy 
which is discussed below [eq. ([25]) for (/x, v) = (x,x)]. 
One can show, however, that the conductivity has no in- 
dependence in the linear order. 

The Hall conductivity a xy in a weak magnetic field is 
obtained in the following way IH 16[ : The Fourier expan- 
sion of the vector potential is chosen as A(r) — A q e lq r . 
Then by the perturbative expansion of eq. (|18[) in terms 
of the Hamiltonian in the magnetic field H—A q -J(—q)/c, 
a xy is given by the three point correlation function, and 
the linear term of the magnetic field B = i(q x A^ — q y A q ) 
is obtained by the q expansion of the Green function 
dk„G = vQa^Q as, 



H 



[q, iv v 



E 

a.=x,y 



A qa 1 

~~chv 



lit, 



At 



an 



dr' 



x e 



' {T T J^q,T)J a {-qy)J v (0^)), 



(25) 



n M1/ (o,u/ TO ) = -i 



B eV 1 



E 



x tr [cr^Q+crxG+cjvGcTyQ — o-^Q + ayQ+a u Qa x Q 
+ (T f _ t G + a x Q + a y Q + a l/ Q — (J f _ l G+o' y G+(T x G+PvG 
+ o^G+OvGoiGayG — &[iG+&vGcryG& x G}- 

for (fJ,,v) — (x, y) in eq. 



(26) 



The matrix trace tr [• 
becomes 

-4[( 1 ^) 2 -( 1 ^+) 2 ][( 1 ^)( 1 ^+ 



(21) [(i£„) 2 - v 2 k 2 



'■[(iuj+) 2 - v 2 k 2 



= F 3 (i£D„,ia)+). 
(27) 
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Using eq. ([22]) , the Hall conductivity is obtained as in the 
same way of the conductivity a: 



3 3 v 2 Bh 



4ir 2 t 



dcj[-f(hj)]A 3 (hw,T,A), (28) 



A 3 (x,T,A) 



' — OO 

.2_ A 2 



4T 2 x\x\ 



7T 

2 ~ 
/(A 2 



-iA 2 



Tan 
-x 2 ) 2 



(T 2 + x 2 + A 2 ) 2 - AA 2 x 2 V 



2rx 



2r|x| 

(29) 



3 

At 



where the last term of A3 which is proportional to 41 
stems from the first and the forth terms in eq. (|22[) 
zero temperature, we have a xy — — ^pr^-Ai (//, T, A). 
The chemical potential dependence of the Hall conduc- 
tivity is shown in Fig. [ljc). In the dirty limit T — > 00, 
and A, /1 — > 0, the Drude-Zener formula is obtained with 
r = H/2t, where r is the mean-free time of qusiparticles. 
This is consistent with the result obtained by formalism 
of the Landau quantization[5]. 

Thermal Hall effect — Finally, we consider the thermal 
conductivity. The thermal transport of the Dirac system 
is discussed in the same way of Refs. 11, 17, 3): For the 
Hamiltonian density 



e(r) = &(r)H V(r), 



(30) 



where the differentiation operator in the one particle 

Hamiltonian (fT]) is replaced as d a — > (d a — d a )/2, by the 
partial integration of eq. |2[). Then the time derivative 
of the local energy is given by 

e = tft# $ + tftfto*- (31) 

The energy current Je is determined so that it satisfies 
the following continuum equation, 



e + V • J E = 0. 



Thus we obtain 

■ J E 



¥<j a y). 



(32) 



(33) 



The current-energy (Tlrtfr ) and the energy-energy (n^f ) 
correlation functions are defined as in the same way of 
eq. (|18p . These correlation functions have similar struc- 
ture of the current-current correlation function: 



lim 



/31m n^f (0,o; + 17?) 



e 2 v 2 B 
Att 2 c 



/3 2 Im Hgf(0,c 
lim =- 

fto; 



dar/'(a:)[/3(x-Ai)]^3(a:,r,A), 

17?) 



(34) 



47T 2 C 



In the zero temperature limit, using the relation 
/'(*)[/% - n)] n dx = 1 (71 = 0), Q (n = odd), 
2n!(l - 2~ ,l+1 )C(n) (n = even > 2), it turns out that 
eq. ([31= and eq. ((35j)= -^f^f ^ 3 (/j, T, A) where 
C(2) = 7r 2 /6 is used. Since the thermopower related 
to does not contribute to the thermal Hall con- 

ductivity n xy , the Wiedemann- Franz law n X y/&xyT — 
(7r 2 /3)(fcB/e) 2 is satisfied for the Hall and the thermal 
Hall conductivities in the low temperature limit. 



Summary — We have discussed the magnetic suscepti- 
bility, the Hall conductivity and the thermal Hall conduc- 
tivity of the 2D Dirac fermions in the weak magnetic field. 
We have shown that strong diamagnetism appears even 
in the gapped system if the Fermi energy is in the middle 
of the gap. We have also shown that the magnetic sus- 
ceptibility derived by the weak magnetic field formalism 
is equivalent with that of the Landau quantization with 
the Euler-Maclaurin formula. The other results includ- 
ing the Hall conductivity are also consistent with those 
of the Landau quantization formalism. 
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